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Abstract 

Bose-Einstein condensation of W bosons in the early universe is studied. It is shown 
that, in the broken phase of the standard electroweak theory, the condensed W bosons form a 
ferromagnetic state with aligned spins. In this case the primeval plasma may be spontaneously 
i— 1 ' magnetized inside macroscopically large domains and form magnetic fields which may be the 

seeds for the observed today galactic and intergalactic fields. However, in a modified theory, 
e.g. in a theory with stronger quartic self interactions of gauge bosons e.g. due to a smaller 
value of the weak mixing angle, antiferromagnetic condensation is possible. In the latter case 
W bosons form scalar condensate with macroscopically large electric charge density i.e. with 
a large average value of the bilinear product of W- vector fields but with microscopically small 
average value of the field itself. 

& . 

i — 1| 1 Introduction 

Bose-Einstein condensation is the quantum phenomenon of the accumulation of identical bosons 

C"\J | in the same state, which is their lowest energy (zero momentum) state. Under these conditions 

they behave as a single macroscopic entity described by a coherent wave function rather than 

p.1 a collection of separate independent particles. Even though the Bose-Einstein condensation 

. — * \ had been foreseen long time ago (1925), it took seventy years to make the first experimental 

observation, which was performed in a dilute gas of rubidium pp. Difficulties in performing this 

observation were created by the extremal conditions necessary for the condensation. Indeed, the 

Bose-Einstein condensation takes place when the inter-particle separation is smaller than their 

de Broglie wavelength, X^b ~ 2ir/\/2mT, so the system must be cooled down to a very low 

^ ■ temperature at ordinary densities. 

In the recent years the study of Bose-Einstein condensate (BEC) became an active area of 
research in different fields of physics from plasma and statistical physics (for a review see book 
[2] and references therein) to astrophysics and cosmology. 

The presence of charged Bose-Einstein condensates has interesting consequences in gauge field 
theories. For instance, in our recent papers [HH], electrodynamics of charged fermions and scalar 
bosons was considered with the fermion asymmetry above a certain critical value such that the 
boson chemical potential had to be equal to the boson mass. In these conditions the bosons 
should condense to ensure electric neutrality of the plasma. The screening of impurities in such 
plasma is essentially different from the case when the condensate is absent. 

A similar problem has been recently considered also in the framework of an effective field 
theory [5j [Q [7j . The authors focused on the applications to the astrophysics of the helium white 
dwarfs discussing possible condensation of helium nuclei and analyzing the thermodynamical 
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properties of the system and possible observational signatures. The formalism they used is 
complementary to ours, but our results agree in overlaping areas. 

In this work we consider the BEC of charged lU-bosons, which may be formed in the early 
universe, and study its magnetic properties. In general W bosons may condense both below 
and above the EW symmetry breaking if the cosmological lepton asymmetry happened to be 
sufficiently high, i.e. if the chemical potential of neutrinos was larger than the W boson mass at 
this temperature. The condensation of vector bosons differs from the theoretically simpler case 
of scalar bosons due to the presence of an additional degree of freedom, their spin states. In both 
cases, scalar and vector, the condensed bosons are in the zero momentum state but in the latter 
case the spins of the individual vector bosons can be either aligned or anti-aligned. These states 
are called respectively ferromagnetic and anti- ferromagnetic ones, see e.g. ref. [2]. The realization 
of one or the other state is determined by the spin-spin interaction between the bosons. In the 
lowest angular momentum state, I = 0, a pair of bosons may have either spin or 2. Depending 
on the sign of the spin-spin coupling, one of those states would be energetically more favorable 
and would be realized at the condensation. In the case of the energetically favorable higher spin 
state, 5 = 2, the vector bosons condense with macroscopically large value of their vector wave 
function (Wj). In the opposite case of the favorable 5 = state the vector bosons form the scalar 
condensate with pairs of vector bosons making a scalar "particle" . Such phenomena are observed 
in solid state physics with such spin-1 bosons as 23 Na, S9 K, and 87 Rb nuclei, see refs. [2], [8j-[ll|. 

Usually experimental studies of the properties of the spin-1 condensate are performed in 
external magnetic fields. Under such conditions the spins of the vector bosons are aligned (frozen) 
due to the interactions of their magnetic moments with the external field. However, in optical 
traps an external magnetic field is absent and spin alignment or de-alignment depends upon 
the internal dynamics of the system. Correspondingly either ferromagnetic or scalar ground 
state would be formed depending on the scattering length of vector bosons in different angular 
momentum channels [8]. 

Let us stress again that in the ferromagnetic case the total spin of the vector bosons has 
macroscopically large value, while in the antiferromagnetic case the lowest energy state of many 
condensed vector bosons has zero total spin, Stot = 0, if their number is even, and has Stot = 1 
(still macroscopically small), if their number is odd. 

Due to macroscopically large value of the total spin in the ferromagnetic spin state the system 
can be accurately described by the mean field approximation, as is argued e.g. in sec. 12.2 of 
book [2] or refs. |X0(, lllj. Indeed, the validity of the mean field approximation is determined by 
the relative magnitude of the fluctuations near the ground state. The fluctuations are induced 
by the particle scattering which can change the spin value in a single reaction by ±1. It is clear 
that for a large value of the total spin the relative fluctuations SS/S ~ 5N/N ~ 1/y/N <C 1, 
while for a small total spin value 5S/S ~ 1. 

In presence of sufficiently strong external magnetic field all the vector boson spins condense 
in the same direction and thus the whole body forms a single magnetic domain independently 
on the spin-spin interactions of the vector bosons, ferromagnetic or not. If however, an external 
magnetic field is absent, several magnetic domains would be formed in the ferromagnetic case, 
due to dynamical instability, and none in the antiferromagnetic case. The discussion of this 
phenomenon in solid state physics and the list of references can be found e.g. in ref. |12j . 

The situation with condensation of lU-bosons is not as complicated as in solid state physics 
where the dominant contribution to the spin-spin interactions comes from the quantum exchange 
effects. However, for the system we are considering here, the exchange forces are not essential and 
the spin-spin interaction is determined by the interaction of the magnetic moments of the vector 



bosons and by their self- interactions. The interactions with relativistic electrons and positrons 
are neglected, see below. 

We show that charged vector bosons would condense in maximum spin states and form clas- 
sical vector field, if only electromagnetic interactions between their spins are taken into account. 
In such a case the spontaneous magnetization at macroscopically large scales would take place. 
On the other hand, the local self-interaction of W creates the spin-spin coupling of the opposite 
sign. In the standard theory the magnitude of this coupling is smaller than that induced by 
the photon exchange, while the exchange of heavy Z-boson does not contribute at all into the 
spin-spin interactions of T^-bosons. Thus the spin-spin coupling is dominated by the interactions 
of the magnetic moments of W. In pure electrodynamics magnetic fields are not screened and so 
one may expect that the plasma effects would not eliminate the dominance of the interactions 
between the magnetic moments. However, the situation is not clear in non-Abelian theories |13j 
and in principle the screening might inhibit the spin-spin magnetic interaction, see sec. 14.41 If 
this is realized, the local quartic selfcoupling of W would dominate over the electromagnetic one 
and W bosons should condense in antiferromagnetic state and form a scalar condensate. Hence 
a classical vector field would not be created. 

The condensation of gauge bosons of weak interactions was considered in the pioneering 
papers [14| . where it was argued that at sufficiently high leptonic chemical potential a classical 
Wj boson field could be created in the early universe. Our results are similar to that of ref. [14] as 
long as the ferromagnetic case is realized. In this case the spins of W add up coherently creating 
classically large average vector Wj boson field. In the hypothetical situation of a stronger quartic 
self-coupling of W we arrive to an opposite conclusion of vanishingly small classical W field but 
with macroscopically large number density of VF-bosons at rest, which is given by the bilinear 
product nw = i(dtWjWj — WWtWj) (as we see in what follows, this expression for the number 
density of W is true in the gauge where the electromagnetic potential is zero). 

A different mechanism of formation of W-boson spin condensate by chaotic magnetic fields, 
which might exist in the early universe, was considered in ref. [15]. If such fields were sufficiently 
strong, this mechanism could operate independently on the spin-spin interaction of VF-bosons and 
would align their spins in the domains with the size of the order of that of the original cosmic 
magnetic field, which are normally microscopically small. To this end a magnetic field with the 
strength B > am^ would be necessary. Such fields might be generated at the electro-weak phase 
transition. The alignment of the W spins reminds the alignment of vector fields in magnetic traps 
mentioned above. Moreover, such an alignment under the influence of a sufficiently strong exter- 
nal magnetic field would take place in both ferromagnetic and antiferromagnetic cases. However 
when the external magnetic field is switched-off or redshifted, the spins of the "antiferromag- 
netic" H^-bosons would be dis-aligned making scalar condensate, while in the ferromagnetic case 
macroscopically large domains with aligned spins would be created. The mechanism of formation 
of such domains is different from the normal ferromagnets due to an absence of the exchange 
forces. So probably the size of the domains is not determined by the usual competition between 
the volume and surface energies but by the causality effects. 

A recent application of vector BEC to astrophysics was studied in ref. [16j . where the con- 
densation of deuterium nuclei was considered. The authors argue that the interaction between 
deuterium nuclei forces them into the lowest spin antiferromagnetic state. 

The paper is organized as following. In sec. 2 kinetic equation approach to BEC is considered. 
In sec. 3 the same phenomenon is studied on the basis of the equations of motion of the vector 
fields. In sec. 4 the spin-spin interaction of W bosons in the zero momentum state is calculated. 
In sec. 5 discussion of the results and conclusion are presented. 



2 Formation of Bose-Einstein condensate. Kinetic approach 

It is convenient to describe formation of BEC using kinetic theory approach. Let us consider 
a system of bosons and fermions in thermal equilibrium with temperature T. As is known, the 
equilibrium distribution functions, up to the spin counting factor, are equal to: 

fF ' B = ex V [(E-VF,B)/T]±l ' (1) 

where the signs plus and minus stand respectively for fermions and bosons and fiF,B are chemical 
potentials of fermions and bosons. Evidently the chemical potential of bosons cannot exceed their 
mass, [is < tub- Otherwise their distribution would not be positive definite. This upper bound 
on hb enforces the Bose-Einstein condensation if the asymmetry between bosons and anti-bosons 
is so high that even maximally large chemical potential, ^b = ttib, is not sufficient to provide for 
such a high asymmetry. In this case the bosonic distribution function takes the form: 

j B = C8^ (p) + 1 , (2) 

exp[(£ - m B )/T\ - 1 

where the first delta-function term describes the condensate and a new constant parameter C is 
its amplitude. The second term, which is the usual Bose distribution, describes non-condensed 
particles and vanishes at T = 0. 

It is easy to verify that the distributions written above annihilate the collision integral, i.e. 
they are the equilibrium distributions. Indeed the kinetic equation has the form: 

^r=Icoii[f], (3) 

where the collision integral is equal to 
T &*)< 

1 COll — nTp 

and 



JdT' m dr fm 5^ (j^Pin -Y,Pfin) \Af\ 2 F[f], (4) 



F[f] = uf in (i ± f fin ) - n/ /in (i ± f in ) . (5) 

In eq. (JU) dr[ n is the phase space of all initial particles except for particle under scrutiny, i.e. 
the initial particle 1, drfi n is the phase space of all final particles, A^j is the amplitude of the 
transition from an initial to a final state. The product is taken over all initial (in) and final (fin) 
states and signs plus or minus stand for bosons and fermions respectively. It is straightforward 
to check that F[f] vanishes on distributions (pQ) for arbitrary fip and [is < tub satisfying the 
usual condition of chemical equilibrium: 

Yl Vin = Yl Win ■ (6) 

If hb = tub, there arises an additional freedom that F[f] vanishes for distribution ([2]) with 
an arbitrary C. The value of the latter is determined by the magnitude of the asymmetry 
between bosons and anti-bosons. Notice that in equilibrium chemical potentials of particles and 
antiparticles are opposite, p, = —fj,. Hence if bosons condense with /i# = tub, the anti-bosons 
cannot condense because \i B = —tub- 

We assume above that the collision amplitude is T-invariant but even if this restriction is lifted 
the collision integral is still annihilated by functions (pQ) or ([2]) due to S-matrix unitarity |17j . 



As is stated in the Introduction, a large lepton asymmetry is necessary to make W bosons 
condense. In particular, the condensate is formed at lepton number densities higher than the 
critical one n c v = m^//67r 2 . When T — > 0, the VF-boson mass, mw, approaches its usual vacuum 
value, created by the Higgs condensate. When the temperature rises, there are two opposite 
effects on inw- The first one is the usual positive temperature correction 5mw ~ eT. The 
second effect is negative and is related to a decrease of the amplitude of the Higgs condensate. 
As a result, at temperatures above the electroweak phase transition, when the Higgs condensate 
disappears [18] . the sole contribution to IF-mass comes from the temperature corrections [19] 
and mw may be much smaller than its vacuum value. 

The temperature dependence of n c v below the EW phase transition was analyzed in refs. [20, 
121] , It should be noted here that, if T is higher than the critical temperature of the EW phase 
transition, the W condensation may take place also for lower, but still large, values of the lepton 
asymmetry [21] . In this case, as is mentioned above, the W mass would be essentially determined 
by the medium effects, i.e. by the temperature corrections or by the condensate itself. Thus the 
W mass in the electroweak symmetric phase might be lower than in the broken phase and a lower 
lepton asymmetry would be required for the condensation. In this paper we mostly assume that 
the temperature is below the electroweak phase transition and thus the plasma consists of massive 
W and Z bosons, neutral Higgs bosons, quarks, leptons, and their antiparticles. Nevertheless the 
interesting possibility of W condensation at lower lepton asymmetry should be kept in mind. 

Models of generation of large lepton asymmetry were considered in refs. |22]-|25], However, 
if the asymmetry was generated above the electroweak (EW) phase transition, it might be trans- 
formed into the baryonic one by the sphaleron processes, creating unacceptably large baryon 
asymmetry So we need to assume that the lepton asymmetry should be created below the EW 
phase transition. 

On the other hand, a mechanism to avoid generation of too large baryon asymmetry could 
be triggered by a large lepton asymmetry itself. It was pointed out in ref. [14] and confirmed 
in several subsequent papers [261 H3] that a large lepton asymmetry suppresses restoration of 
the electroweak symmetry and hence if lepton asymmetry was generated at very high tempera- 
tures when sphalerons were not active and the electroweak symmetry was broken at these high 
temperatures, the electroweak baryon non-conservation would never be efficient. 

Finally, concerning the formation of the Bose-Einstein condensate, it is essential that the 
particles in question possess a conserved quantum number that forces their chemical potential 
to be nonvanishing, if the number of particles is not equal to the number of antiparticles. The 
amplitude C is then calculated from the known difference of the number densities of particles 
and antiparticles. 

In this paper we consider a simple example of electrically neutral plasma in which charged W- 
bosons condense because of a large asymmetry between leptons and antileptons. For simplicity 
we confine ourselves to only one family of leptons. This simplification does not influence the 
essential features of the result. A more detailed analysis with all the particles included can be 
found in ref. [27]. Quarks may be essential for the imposing of the condition of vanishing of all 
gauge charge densities in plasma and for the related cancellation of the axial anomaly but we 
work in the lowest order of the perturbation theory where the anomaly is absent. 

The plasma is supposed to be electrically neutral, with zero baryonic number density but 
with a high leptonic one. The essential reactions are the direct and inverse decays of W: 

W + o e + + v . (7) 

The equilibrium with respect to these processes imposes the equality between the chemical po- 



tentials: 

HW — \iv — He (8) 

The condition of electroneutrality reads: 

n w + - n w - - n e - + n e + = (9) 

The leptonic number density is equal to 

til = n v — n p + n e - — n e + . (10) 

Here n = g s J d 3 p/(2ir) i f is the number density of the corresponding particles and g s is the spin 
counting factor. One should remember that only left-handed electrons participate in reaction ((7|), 
chirality is conserved in reactions with Z-bosons and photons, and chirality flip may take place 
only in reactions with Higgs bosons. 

Leptonic number is supposed to be conserved, so ul is constant in the comoving volume and 
is a fixed parameter of the scenario. (B+L)-nonconservation induced by sphalerons is neglected 
because temperature is mostly assumed to be below the electroweak phase transition. 

It is evident that for sufficiently high hl the chemical potential of W should reach its maximum 
value, fj>w = mw-i an d with further rising ni, W-bosons must condense. 

Let us first assume that the density of leptonic charge is very large, n^ > T 3 . Correspondingly 
the amplitude of the condensate must also be large, n w « C/(2ir) 3 > T 3 . In this case \x e ^ v > 
mw and equations (fS HTUj) can be easily solved: \i v ss ^g, m ss ^ 3 /2-7r 2 , and n w w (2/3)riL- 
Interparticle separation of W-bosons under these conditions is d ~ n L < T _1 , while the de 
Broglie wave length of the condensed W's with zero momentum is formally infinitely large. In the 
realistic condensate the particle momentum is not precisely zero but is of the order of the inverse 
size, d , of the region where IV-bosons condense. Another relevant quantity is the de Broglie 
wave length of W-bosons above the condensate, A^^ ~ 2ir/\ / 2mfT. The condition Xjb > d would 
be fulfilled as long as T < my/ /2tt 2 . 

A huge cosmological lepton asymmetry, ul 3> T 3 , could be created in the Affleck-Dine sce- 
nario [28J. In this model the universe could be quite cold. Later when Ul ~ 1/a 3 is diluted by the 
cosmological expansion down to the value when fiyy becomes smaller than mw the condensate 
would evaporate and the universe would be reheated. 

Though the possibility of a huge lepton asymmetry is quite interesting, the condensation of 
W-bosons could take place even with much smaller ul ~ m w-> as can De seen from the analysis 
of the kinetics presented in this section. The interparticle separation of W-bosons under these 
conditions is d ~ n L ~ m w . In this case the temperature may be relatively high, T ~ my/, but 
nevertheless VF-bosons could condense. According to the equilibrium distribution (|2|) the plasma 
would consist of two parts, the condensate with zero momentum and the high temperature plasma 
over the condensate. The de Broglie wavelength of the high temperature plasma would be again 
A^ ~ 2-7r/v / 2mT, that is larger than d as long as T < 2-K 2 my/- As is argued above, the de Broglie 
wave length of the condensed W-bosons even in this case is much larger than the interparticle 
distance. 

Thus at T ~ my/ the charged weak bosons, W, condense if all the relevant quantities are close 
by an order of magnitude to the VF-boson mass to a proper power, i.e. [i v ~ \x e ~ raw-, ni ~ m w . 
They rise with rising temperature, as can be found from numerical solution of eqs. (J81 110I) . 

We will also discuss the W condensation above the EW phase transition. In this situation the 
Higgs condensate is absent and the M^-bosons acquire, due to temperature corrections, mw ~ 
gT |19| . where g ~ 0.1 is the typical coupling constant. 



The critical values of chemical potential necessary to achieve W condensation at lower m, 
have been calculated, even though an analytic solution is possible only in some special limits 
and the approximations used are not valid in some regions of parameters. The general equations 
together with some critical discussion concerning the approximations used can be found in the 
literature, see |2CH 121 j . For instance, at high temperature, that is T much larger than masses 
and chemical potentials but smaller than the EW breaking value, VF-bosons would condense if 
Hl > 0.3\/T 2 — T 2 , where //£ is the chemical potential associated to the total lepton number and 
T c ~ 200 GeV [21] . Another limit that can be analytically calculated is the low temperature one 
(T ~ 0). Under this hypothesis, VF-bosons condense when //£ > 88.4 GeV +T 2 /64.8 GeV. 

A caveat for a large lepton asymmetry arises from the consideration of the big bang nucle- 
osynthesis with strongly mixed neutrinos. It is shown [29] that leptonic chemical potentials of 
all neutrino flavors are restricted by \fi u /T\ < 0.07 at the BBN epoch. However, it should be 
noted that the entropy release from the electroweak epoch down to the BBN epoch diminishes 
the lepton asymmetry, n^/n-y, by the ratio of the particle species present in the cosmological 
plasma at these two epochs, which is approximately 10. Hence the original lepton asymmetry, 
{ni/n^EW could be of order unity and still compatible with BBN. It should be also noted that 
the BBN bound can be avoided [30] if neutrinos are coupled to the hypothetical pseudogoldstone 
boson, Major on. 

Moreover, at high temperatures, when the Higgs condensate is underdeveloped, the W-hosoxi 
mass may be noticeably smaller than the plasma temperature and W could condense even at 
|/i/T| < 0.07. If Ws condensed above the electroweak phase transition, the magnitude of the 
chemical potential necessary for the condensation was much smaller than temperature, nw ~ gT, 
as explained above. As we see, in this case the chemical potentials of electrons and neutrinos 
would be also much smaller than temperature. In this limit the differences between number 
densities of fermions and antifermions are: 

n F -n F = -g F /i F T 2 , (11) 

D 

where g F is the number of the spin degrees of freedom, g e = 2, g v = 1. 

Substituting this expression into eqs. ([HV llOp and assuming an arbitrary chemical potential 
H\y < T, we find that the condensate would be formed if approximately til > gT 3 and that all 
the chemical potentials are of the order of gT. Correspondingly the necessary lepton asymmetry 
could be as small as til/ti^ ~ g. We see that even without the entropy release such lepton 
asymmetry is not dangerous for BBN. 

3 Equations of motion of gauge bosons and their condensation. 



In this section we consider the evolution of the gauge bosons of SU(2) x U(l) electroweak model 
and analyze the possibility of the condensate formation as supplementary to the approach of the 
previous section. In addition to that, using the equations of motion allows to take into account 
the spin-spin interactions of the vector bosons and to check if the condensate is ferromagnetic or 
antifer r omagnetic . 

The Lagrangian of the minimal electroweak model has the well known form: 

L = L g b + L sp + L sc + Ly u k, (12) 



which are respectively the gauge boson, the spinor, the scalar, and the Yukawa contributions. 
Explicitly we have: 



1 ■ 1 

_/~n fti fJ.u _ f f fiu 



(13) 



^"av ~ duA v o u A + ge A A , 



fi vi 



fuu = duB u - d v B i 
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i = 1,2,3 



L 



sp 



mpw, 



d^ 



d ^- l -gaiA^- l -g'YB,)AK 



L sc = - (D^ (D M $) + -» 2 &$ - -\(&*) 2 , 



D^ 



d^-'-g^A^-'-g'B,)^. 



and the Yukawa Lagrangian describes interactions of fermions with the Higgs field. 

In the expressions above A l „ and B^ are the gauge boson potentials of the SU(2) and £7(1) 
groups respectively, g and g' are their gauge coupling constants, Y is the hypercharge operator 
corresponding to the U(l) group and a 3 are the Pauli matrices operating in SU(2) space. As 
usually, the repeated indices imply summation. 

In the broken phase the physical massive gauge boson fields are obtained through the Weinberg 
rotation: 



W, 



± 



V2 



c wA u - swBn , 



Au = s w Al + c w B : 



A" 



(14) 



where cw and sw stand respectively for cos 6w and sin 6w and 6w is the Weinberg angle. The 
other fields involved in the Lagrangians presented above are the spinor \& and the Higgs field 
$ = [(/> + ,(/>°] . The latter, after the symmetry breaking, acquires non-zero vacuum expectation 
value, v, and takes the form 3> = (1/V2)[0, v + (p ( [\ T , where the upper index T means transposition 
and 4>i describes excitations in the broken symmetry phase, i.e. neutral massive Higgs particle. 
We are considering here one lepton family but the results can be easily generalized to the three 
family model. We use the unitary gauge in which the particle content is explicit, for example 
physical gauge bosons have three polarization states and only one physical neutral Higgs field, 
4>\ is present. 

From the presented above equations one can conclude that in addition to the usual kinetic 
and mass terms, vector boson interactions contain the following cubic and quartic couplings, see 
e.g. ref. [E]: 



L 3 = iecot W [(d^W v - d v W») W\Z V - (d^W* - d v W^J W^Z U + W^Wl {d»Z v - d u Z») 
+ ie Ud^Vr - d u W fl ) W},A V - (Vir 4 - d v W^) W^A U + W^Wl (d»A v - d v A")] ,(15) 

2 r 

% — (wlwrf - w\w^w v w 

2sm 2 6 w [\ » J » 

- e 2 cot 9 W \2W\W^Z V A V - W ] U Z^W V A V - W\A*W V Z 



L 



e cot 9\y 



wlw^ZvZ" - wlz»w v z v 



WlW»A v A v - W\A^W V A V 



(16) 



Using the standard Euler-Lagrange procedure we can obtain the following Maxwell equations 
for the electromagnetic field: 



d„F\ 



Jf, + ie \w},d u W^ - W d v Wl + Wl d u W^ - W v d u W^ + 2W"d a Wl - 2W^d u W l 



v "p- 



v '-'p ' 



+ 
+ 



2A v WlW» - A^ (wlW v + W^Wl 



e cot 9w 



2Z u wlw> i - z" (wpr v + w^wl 



(17) 



where A„ is the electromagnetic potential, F^v = duA v — d u A^ is the Maxwell tensor, and jf is 
the electromagnetic current of the charged fermions. All the rest in the r.h.s. of eq. (|17p can be 
understood as the electromagnetic current of W bosons, J^J . In what follows we assume that the 
total electric charge density of the plasma is zero and the average three-vector current vanishes 
as well. 

Equations of motion for the other vector fields , W^ and Z, can be obtained from the 
Lagrangians (|15l [16]) plus the contributions from the kinetic and mass terms: 



d^W? + mivWt 



W yy v 



= ie [A»W^ U - d^WAv) + d^WuA") - W»F^ V ] 

+ ie cot 9 W [Z^W,,,, + dp(Z»W u ) - M (W"Z v ) + WZ VV \ 



+ 



(eV 



8^ + m\Z v 



sin 2 W ) [Wi,(WjW M ) - Wl{W^W^)\ + e 2 cot 2 9 W {W V Z^ - Z V Z„W») 
+ e 2 cot 9 W (2W U Z^ - Z U W^ - A U W^Z^) + e 2 {W v A»A„ - A U W^), (18) 
= ie cot e w [w^wl„ - W^W^ + d^W^Wl) - d^W^W v ) 

+ e 2 cot 2 6 W {2Z u WlW^ - WlWpW* - W^W^Z^ 

+ e 2 coX6 w (2A v wlW^-W} l W ll A^-W v WlA l )j, (19) 

where V^ = d^V u — d v V^, V^ = W M , Z^ and we assumed that the fermionic currents coupled to 
W and Z are absent. To avoid confusion it is worth noting that the field strength defined after 
eq. (fT3|) includes the gauge boson self-coupling, while it is not included into V^. 

Equations (118|19p are the equations for the corresponding field operators. They describe 
classical fields in the tree approximation and do not include the effects of W and Z instability. 
The latter can be taken into account by perturbative iteration of these equations. The effects of 
instability are discussed below. 

We will consider the case when the electric charge density of fermions, Jq , is non- vanishing 
and homogeneous. It is usually assumed that the primeval plasma is electrically neutral and thus 
the non-zero charge density of fermions must be compensated by the opposite sign charge density 
of W. To study such a system it is convenient to use the electromagnetic gauge freedom and to 
impose the condition A^ = 0. We also assume that the average value of Z^ = 0. In this case 
there exists a homogeneous solution of the equation of motion of the form: 



Wf, = C^exp(-im w t), 



(20) 



where we impose the condition d^W^ = to eliminate the non-physical spin degrees of freedom. 



So Ca is a constant vector with vanishing time component and m 
by the nonzero vacuum expectation value of the Higgs field 



nw is the boson mass determined 
Q In contrast, 



m some papers 



4 We want to stress an important difference between Bose-Einstein condensate and Higgs condensate, which is 
often overlooked. The equation of state of the former is simply P = 0, while for the latter P — —p, where P and 
p are respectively pressure and energy densities. 



the gauge condition of time independent charged vector field is taken: W^ = C^, while the 
electromagnetic vector potential is non-zero: A^ = fi^mw/e- 

In addition to the Higgs induced part, the mass of Wn contains also the contributions induced 
by the temperature effects [19] and by the impact of the VF-condensate itself, which are disre- 
garded at this stage, see discussion in sec. [2] about possible condensation of W at lower lepton 
asymmetry. 

Solution ([20]) corresponds to the Bose condensate of VF-bosons describing a collection of these 
particles at rest. This field configuration corresponds to the condensation of the positively charged 
W + . The condensation of W~ is described by the complex conjugate expression. However, such 
a solution is not obligatory and depends upon the kinetics of the system and the interactions 
between IF-bosons at rest. In fact the only condition which we have at this level is the condition 
of the electric neutrality and it demands only that the average values of bilinear combinations 
of W must be non-zero, while classical vector field Wu may vanish on the average. A possible 
vanishing of the classical vector field Wj, where j = 1, 2, 3 is the spatial vector index, is physically 
quite evident. The condensate is a collection of VF-bosons at rest with the charge density which 
compensates the charge density of fermions. The most probable state of such particles (the highest 
entropy state), if the spin-spin interaction is neglected, is the state with chaotic distribution of 
the individual spins. It is natural to expect that the average value of the total spin in such a state 
is zero or at least not macroscopically large. The situation is opposite in the ferromagnetic case 
when the spin alignment is energetically favorable and the classical vector field could be formed. 

In reference [14J a similar statement of the formation of a classical vector field was done but 
without an analysis of the dynamics introduced by the spin-spin interaction. In the quoted paper 
the mentioned above gauge of time independent W is used, but the physical results are of course 
gauge invariant. In this gauge the condition of the charge neutrality becomes 2AqW^W~ = —Jq, 
which is again bilinear in W field and from the condition of non-zero charge density of the 
condensed W's does not follow that there exists the classical field Wj ^ 0. 

One more comment may be in order here. The instability of W can be introduced in the usual 
way by adding an imaginary part to the mass equal to the decay half-width. The introduction 
of such a term into the equation of motion for W leads to the exponential decay of the field, 
W ~ exp(— Tt/2). However this description is valid only for the decay into vacuum. 

For the decay into a dense medium the Fermi exclusion principle should be taken into account. 
Hence, if neutrinos have sufficiently large chemical potential, such that all the states where W 
could decay would be occupied, the decay rate would be exponentially suppressed and solution 
(|20l) could be physically realized. This observation establishes the equivalence between the kinetic 
approach of sec. [2] and that presented here. 

In fact the absolute stability of the condensed W s is unnecessary. Even if IF-bosons decay, the 
ferromagnetic condensate can be formed, if the time of the condensate formation is shorter than 
the life-time of VF-bosons in the plasma. The condensed VF-bosons are in the state of dynamical 
equilibrium: VF's evaporating from the condensate because of their decay or scattering of hot 
fermions, are compensated by VF's coming back by the inverse processes. In thermodynamical 
equilibrium the average number of VF-bosons in the zero momentum state remains constant. The 
decay rate of the VF-bosons in plasma is proportional to 

r ° (21) 



[ e -(«W2+Me)/T + ]_] r e -(m w /2-n„)/T + 1] ' 

where Tq ~ amw is the decay rate of W in vacuum. The denominator in eq. (|2ip comes from 
the Fermi suppression terms (1 — / e +)(l — f v ). We consider the decay W + — > ue + and take 
into account that in thermal equilibrium the chemical potentials of electrons and positrons are 
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equal by magnitude but have the opposite signs. In the case of very large lepton asymmetry, 
nL, H>v 3> T, the decay rate of VF-bosons would be exponentially suppressed and their life time 
can be longer than the time of the spin alignment, t$- As we mentioned above, the cosmological 
generation of L 3> T 3 may be realized in a version of the Affleck-Dine scenario [28] which leads 
to a cold universe with non-relativistic PF-bosons. Moreover, even in absence of the exponential 
suppression and relatively small lepton asymmetry, L ~ T 3 , the life-time of VF-bosons in the 
plasma can be larger than rg. As is shown below, the Hamiltonian of spin-spin interaction is 
given by eqs. (I29|34p . Correspondingly the characteristic time of the spin alignment can be 
estimated as: 

T S ~ 1/U s ? m ~ m 2 w /(n w e 2 S 2 ) , (22) 

where S is the total spin of the condensed W-bosons, ny/ = 1/d 3 is their number density, and 
d is the average distance between them. We approximated <5(r) as l/<i 3 = ny/. Evidently ts 
can be considerably shorter than the life-time of VF-bosons. The same "stability" arguments are 
applicable to the decay of W into quarks. 

4 Spin-spin interactions of W bosons 

As we mentioned above, the form of the vector boson condensate depends upon the interaction 
between the vector bosons at rest. If the latter favors the opposite spin configuration, i.e. a pair 
of bosons "prefers" to be in the zero spin state, the condensate would have zero total spin, i.e. 
PF-bosons would form the scalar condensate (antiferromagnetic case). In the opposite situation 
of the favorable spin-two state, the spins of all vector bosons in the condensate would be aligned 
and the condensate would have macroscopically large spin (ferromagnetic case). 

First, in subsection 14. II we will consider only the spin-spin interaction induced by the electro- 
magnetic interactions of VF-bosons, namely by the coupling of their magnetic moments. Next, 
in subsection 14.21 we will take into account the local quartic self-coupling of W. In subsection 
14.31 we take into account the spin-spin interaction caused by the Z-boson exchange between Ws. 
In sec. 14.41 we discuss the effects of possible screening of the magnetic moment interactions by 
plasma effects. In section [431 we discuss the results of the previous subsections and conclude. 

4.1 Electromagnetic interactions 

The essential particles in the system we study in this section are the charged W bosons at rest and 
charged relativistic fermions, which ensure the electric neutrality of the medium. The latter are 
electrons (or positrons) and quarks but these details are not important. For relativistic fermions 
helicity is conserved and hence the interaction of their spins with the spins of W is not essential, 
because on the average the electron-positron medium is not spin-polarized. Accordingly we take 
into account only the spin-spin interaction of non-relativistic VF-bosons and disregard the impact 
of the charged fermions. The electromagnetic interaction between VF-bosons is similar to the 
well known interaction of nonrelativistic electrons, which is described by the Breit equation. The 
derivation of this equation can be found e.g. in book |32j. According to the Breit equation the 
interaction between the magnetic moments of two electrons (i.e. their spin-spin interaction) has 
the form: 

U"lr) - J^ [<^> - "(■* ■'>(■»■" - * (ol . ^<3>( r) l , (23 ) 
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where 01 2 are the spin operators of the electrons, i.e. the Pauli matrices averaged respectively 
over the first and second electron wave functions, and e 2 = 47ra = 47r/137 (to avoid confusion let 
us note that in ref. |32j the notation is different, namely e 2 = a). 

The expression for potential (|23p created by the one photon exchange is true for virtual 
photons propagating in vacuum. The presence of plasma changes the propagator and could 
modify the spin-spin potential. This is discussed in section I4.4L The screening effects are not 
important, at least in some temperature range. 

The analogue of the Breit equation for VF-bosons can be derived along exactly the same lines 
as is done for electrons. The electromagnetic interaction between two W bosons in the lowest 
order in the electric charge, e, is described by the usual one-photon exchange diagram. In the 
Feynman gauge, where the photon propagator is D^ u = g^ v /q 2 , the amplitude corresponding to 
this diagram is: 

M = — L^w^w^v ( puP2 , q ) V »( P3 ,p 4 ,-q)W a WP (24) 

{P1-P2Y 

where V a pu is the most general CP invariant W^W^f vertex |33| : 

V a/ 3 fl (pi,P2,q) = ie 



P(i9afi + 2 {qpgan ~ <?a3/3 M ) + ( l ~ M (^5a/i - <Ja90n) + ^Hf" V\A*<lp 



2"% 



,(25) 



where p\ and p% are the momenta of the incoming particles, P2 and p 4 are the momenta of the 
outgoing particles and p = pi + P2, Q = P2 — Pi- This expression should be symmetrized with 
respect to the interchange of VF-bosons in the initial and/or in the final states. 

The vertex written above contains two anomalous coupling parameters fey and A 7 . As we 
can see from eq. (fl7| ). the standard electroweak model predicts, up to the second order in the 
electromagnetic coupling constant e: fey = 1, A 7 = 0. In what follows we mostly assume that 
these values are true, since they are compatible with the present experimental data for triple 
gauge boson couplings [HI]. In this case the amplitude (041 reduces to: 



M = ^W^WP [ Pll g aa , + 2(q a g a/fl - q al g a „)] [jfg m - 2{q m , " - q pi g ")] W?W$ . (26) 



The spin-spin interaction is contained in the product of the last two terms in the square 
brackets in eq. (j26|) i.e. the terms containing vector q. The spin operator of vector particles is 
defined as the generator of the rotation group belonging to its adjoint representation and is equal 
to the vector product: 

Si = -i Wj, x Wi (27) 

Hence the scattering amplitude induced by the interaction between the magnetic moments of the 
charged vector bosons is equal to: 



2„2 



e p 



M s = f^ t (Si • S 2 ) - (Si • q) (S 2 ■ q) (28) 

where p is the ratio of the real magnetic moment of W to its value predicted by the standard 
electroweak theory (e 2 /m 2 v ) and we divided by 4m 2 y for proper normalization of the W-w&ve 
function, as is explained below, see Sec. 14.21 
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The potential which describes the electromagnetic spin-spin interaction is the Fourier trans- 
form of amplitude (f28l) and is equal to: 

«"('•> - ^r [^ " 3 (Sl -^ S2 - r) - f (S, • S 2 ) { «(r)] . ,20) 

This potential has the same form as the corresponding one in the Breit's equation for electrons 
(|23p but with the different numerical coefficient. 

To calculate the contribution of this potential into the energy of two VF-bosons at rest we 
have to average it over their wave function. In particular, in the condensate case, it is a S-wave 
function that is angle independent. Hence the contributions of the first two terms in eq. (|29p 
mutually cancel out and only the third one remains, which has negative coefficient. Thus the 
energy shift induced by the spin-spin interaction is equal to: 

5E= [*L umn {r) = -llfL (Sl • S 2 ) , (30) 

J V 3 Vmyy 

where V is the normalization volume. 

Since S 2 ot = (S± + S2) 2 = 4 + 2S1S2, the average value of S1S2 is equal to 

S 1 S 2 = Sl t /2-2. (31) 

For Stot = 2 this term is S1S2 = 1 > 0, while for Stot = it is S1S2 = — 2 < 0. Thus, if 
the spin-spin interaction is dominated by the interactions between the magnetic moments of W 
bosons, the state with their maximum total spin is more favorable energetically and VF-bosons 
should condense in the ferromagnetic state. This could lead to the spontaneous magnetization in 
the early universe. 

4.2 Quartic self-coupling of W 

The contribution to the spin-spin interactions of W comes from the first term in Lagrangian (1160 
or from the third term in the r.h.s. of equation of motion (|18p . The first term in Lagrangian (|16p 
can be rewritten as: 

\wlW^) 2 - W\W^W V W V \ = % fw f x w)\ (32) 

J 2 sin 6y/ v ' 



Liw - 



2 sin 6w 

It is assumed here that d^W^ = and thus only the spatial 3-vector W is non- vanishing, while 
W t = 0. 

Since the corresponding Hamiltonian is obtained from L^w by changing sign and since spin 
operator (|27p contains imaginary unit factor, the sign of the Hamiltonian is positive. It means 
that the low spin states are energetically favorable. 

For the comparison of this Hamiltonian with potential energy (|29h we need to properly nor- 
malize the wave functions of W. In the Hamiltonian the usual relativistic normalization is used, 
according to which the number density of W is equal to nw = 2my^W' W, while in the non- 
relativistic Schroedinger equation the wave function is normalized to unity, 

d 3 r#| 2 = 1 (33) 

Accordingly the Hamiltonian should be divided by 4m^ and its Fourier transform producing the 
spin-spin interaction potential would be: 

U^ = R 2 e \» (SiS 2 )^)(r). (34) 

arriyy sim &w 
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Thus the quartic self-coupling of W contributes only to the spin-spin interaction whose sign is 
antifer r omagnetic . 

The same result can be obtained from equation of motion (|18p if one takes into account that 
in the non-relativistic limit: 

8 2 W = (-E 2 + m 2 w )W « 2miyeW, (35) 

where e = (E — mw) is the non-relativistic energy. 

4.3 Z-boson exchange 

The contribution to the spin-spin potential between a pair of W from the Z-boson exchange can 
be found from eq. (|18p where we substitute the expression for Z u taken from eq. (|19p in the limit 
of vanishing four- momentum of Z. Indeed the transferred momentum is much smaller than mz, 
and so the diagram with Z-exchange is effectively local with Z-boson propagator equal to \/m 2 z . 
Hence the contribution from the Z-exchange in the e 2 order to eq. (|18p is: 

dpWt + m 2 v W l + 4e 2 cot 2 6 w (mw/mz) 2 (W^W)Wi + ... = , (36) 

where j = 1,2,3 is the spatial vector index. 

We see that the Z-boson exchange does not contribute to the spin-spin interactions of W. 
However, it should be kept in mind that this result is true only for the non-relativistic Z-bosons, 
while above the phase transition the contributions of Z bosons and photons are similar. 

4.4 Plasma screening 

In plasma the time-time component of the photon propagator is modified as 1/q 2 — > l/[q 2 + 
noo(q)], where noo(q) is the photon polarization operator in plasma. Usually noo = m D , where 
mo is the Debye screening mass, which is independent on q. So the pole at q 2 = shifts to an 
imaginary q leading to the well known effects of the Debye screening. As it was found recently [3] 
- [TJ, the presence of the charged Bose condensate drastically changes the polarization operator 
leading to an explicit dependence of noo on q which gives rise to infrared singular terms. The 
modification of the propagator takes place already in the lowest order in the electromagnetic 
coupling, e 2 , i.e. in the one loop approximation. 

On the other hand, the space-space component of the propagator remains massless, Hij ~ 
1/q 2 . It is known to be true in pure Abelian electrodynamics in any order of perturbation 
theory, while in non- Abelian theories the screening may occur in higher orders of perturbation 
theory due to the infrared singularities [13]. The screening may potentially change the relative 
strength of the electromagnetic spin-spin coupling, which is affected by screening effects, with 
respect to local W A coupling which is not screened. However, in the broken phase the system is 
reduced to the usual electrodynamics, where screening is absent and T^-bosons would condense 
in the ferromagnetic state. In the unbroken phase of the electroweak theory the answer is not 
yet known. Perturbative calculations are impossible because of the violent infrared singularities. 
Maybe lattice calculations analogous to those done in QCD will help. 

The potential describing the magnetic spin-spin interaction is related to amplitude (|28l) with 
a modified photon propagator. Thus it can be written as: 

U^Hr) = -"X I i^T^P^L W (Si • S 2 ) - (Si • q) (S 2 • q)l , (37) 



m 2 w J (2vr)3(g2 + n ss (q)) 
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where II SS is the plasma correction to the space-space component of the photon propagator. 

If, as above, we assume that the wave function of VK-bosons is space independent and average 
this potential over space, we obtain the following expression for the spin-spin part of the energy 
shift: 

Clearly SE vanishes if Ii ss is non-zero at q = 0. Of course, this is an unphysical conclusion, 
because the integration over r should be done with some finite upper limit, r max = I, presumably 
equal to the average distance between the W bosons. So instead of the delta-function, <5( 3 )(q), 
we obtain: 

rl 47T 

/ d rexp(iqr) = —^ [siniql) — qlcosiql)} . (39) 

Jo q 

The energy shift is given by the expression: 

e V QiQ j f d 3 q [sm{ql) - ql cos{ql)} [q 2 5ij - q l q j ] 



where V = 4irl 3 /3. 

When we average over an angle independent wave function, e.g. S-wave for the condensate, 
the non-vanishing part of the integral in Eq. (140 h is proportional to the Kronecker delta, hence: 



5E = S{S J 2 A5 ij , (41) 

where the coefficient A can be calculated by taking trace of Eq. (I40p : 



Hence the energy shift of a pair of W-bosons in S-wave state due to the spin-spin interaction is: 
ap fo a ^ 8ne2 P 2 f d3q t sin ( g/ ) ~ ql cos (gQ] fA o\ 

6E -- {Sl - S2) w^Jwr qtf+n s M] ' (43) 

Introducing the new integration variable x = ql we can rewrite it as: 

4g2„2 /-oo ^j, _ -, 

5E = -(Si -S 2 ) 2 / 2 , ;2tt i n\ x sinx + / 2 n sg (x//)cosx , (44) 

^TrVm^y Jo x z + l z n ss (x/l) L J 

We used here the usual regularization of divergent integrals: ex.p(±iql) — > exp(±iql — eq) with 
e — t- 0. With such regularization /^° dxcos(x) = 0. 

Evidently, if II SS = 0, we obtain the same result as that found in sec. 14.11 In fact the nec- 
essary condition for obtaining the "unscreened" result is Z 2 n ss (x/7) <^ 1, but for a large l 2 Ii S s 
the electromagnetic part of the spin-spin interaction can be suppressed enough to change the 
ferromagnetic behaviour into the antiferromagnetic one. This might take place at high tempera- 
tures above the EW phase transition when the Higgs condensate is destroyed and the masses of 
W and Z appear only as a result of temperature and density corrections and thus are relatively 
small. The quantitative statement depends upon the modification of the space-space part of the 
photon propagator in presence of the Bose condensate of charged W. As far as we know, this 
modification has not yet been found. 
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4.5 Discussion 

It is instructive to compare the magnetic properties of solid state bodies with the considered here 
(anti-)ferromagnetism of BEC of VF-bosons. Magnetic properties of matter are determined by 
the state of outer (unpaired) electrons. A pair of electrons belonging to different atoms may be 
either in symmetric, Stot = 1 ; or antisymmetric, Stot = 0; s pi n state. Since the total wave function 
of two electrons must be antisymmetric, their spin state has opposite symmetry with respect to 
their orbital wave function. Symmetric and antisymmetric electron states evidently have different 
energies, which we denote as E s and E a respectively. Accordingly the spin Hamiltonian can be 
written as 

H spin = -JS 1 -S 2 . (45) 

The quantity J = (E s — E a ) is usually called the exchange energy. Its sign is determined by the 
atomic ground state structure. Evidently J > favors parallel spins, while J < favors antipar- 
allel spins. In atomic systems the exchange energy at small distances is typically of the order of 
fractions of eV, that is about 10 3 times larger than the typical direct magnetic dipole (spin-spin) 
interaction between electrons. Hence the exchange interaction may force the magnetic dipoles of 
electrons to be aligned or anti-aligned independently on their direct magnetic interaction. The 
situation changes at large distances, since the exchange energy decays exponentially, while the 
magnetic interaction behaves like r~ 3 . Thus the latter dominates on macroscopic scales, lead- 
ing to the breaking of the system into domains with different directions of the magnetic field 
and consequently to zero net macroscopic magnetization. Nevertheless, the ferromagnetic nature 
emerges when an external magnetic field is applied to the system. 

Fortunately, the system we are considering here is much simpler. All W-bosons are in the 
same state with zero momentum and are not binded into a complicated atomic system. Evidently 
the exchange forces are not essential in this case. Indeed, W^-bosons are in symmetric orbital 
state. Hence their spin wave function should also be symmetric and both the allowed spin states 
of W, the scalar, S = 0, and the tensor, S = 2, ones are also symmetric. 

In principle, electrons and positrons could distort the spin-spin interactions of W by their spin 
or orbital motion and thus destroy the attraction of parallel spins of W. However, it looks hardly 
possible because electrons are predominantly ultra-relativistic and they cannot be attached to 
any single W boson to counterweight its spin. The low energy electrons cannot be long in such 
a state because of fast energy exchange with the energetic electrons. The scattering of electrons 
(and quarks) on VF-bosons may lead to the spin flip of the latter, but in thermal equilibrium this 
process does not change the average value of the spin of the condensate. 

5 Conclusion 

We have studied here the Bose-Einstein condensation of charged weak bosons, W^ 1 . Such con- 
densation might occur in the early universe if the cosmological lepton asymmetry was sufficiently 
large. The primeval plasma is supposed to be electrically neutral and to have zero density of 
weak hyper charge. 

In general a Bose-Einstein condensate of vector fields may form either a scalar state, when 
the average value of vector Wj is microscopically small, or a classical vector state when all the 
individual spins of the condensed vector bosons are aligned at a macroscopically large scale. 

In solid state physics the realization of one or the other of these two possibilities depends 
on the short range exchange forces which dominate over the direct spin-spin interaction at small 
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distances. In the case of W bosons the choice between ferromagnetic or antiferromagnetic state 
is determined by the spin-spin interaction of the individual VF-bosons, realized through the in- 
teractions of their magnetic moments and their quartic self-coupling. 

The total spin-spin interaction potential for W is the sum of two terms (|29l) and (|34l) . As we 
have discussed below eq. (|29p . the first two terms in the interaction of the magnetic moments 
cancel each other after averaging over a S-wave function. Thus only the <5-function term survives. 
In the standard electroweak model p = 1 and thus the absolute value of the coefficient in front 
of SiS2<5' 3 ^(r) in eq. (129p is 2e 2 /(3m 2 y) which is larger than the corresponding term in eq. 
(|34p . Hence the former dominates and the energetically favored configuration of a multi-W state 
should have a macroscopically large total spin. However, as we have pointed out in Sec l4.41 the 
plasma screening of the interaction between the magnetic moments may be dangerous for the 
W-ferromagnetism. In the broken phase the problem is reduced to that of pure QED, where 
it is known that magnetic forces are not screened. However, in non-Abelian gauge theories the 
absence of screening is known only in the lowest order in perturbation theory, while higher order 
calculations suffer from strong infrared divergences and are not reliable. For the resolution of this 
problem non-perturbative methods, as e.g. lattice calculations, are necessary. At the moment 
the problem remains unsolved. 

It should be also noted that, even though Eqs. f|23|) and (|29p are calculated for electrons 
and W bosons, they are valid respectively for any spin-(l/2) and spin-1 species having the usual 
electromagnetic interactions. So they may be applied to other particles, including the ones present 
in the extensions of the standard model. If these particles are in a S-wave state, on the average 
the only delta term survives and hence they may form ferromagnetic states. Indeed, as one can 
see from Eq. (13ip . the lowest energy state for a S-wave function is the one with the maximum 

c2 

D tot- 
It is clear from Eqs. (|29p and (|34p that the condensate of H^-bosons would be antiferromag- 
netic if the bosons have a negative non-standard contribution into the magnetic moment, such 
that p < 3/(16 sin 2 6\y). The ferromagnetism of W can be destroyed also in a model with a 
smaller value of the Weinberg angle. All that demands a strong deviation from the standard 
model and most probably is excluded, but these effects may be important in applications to 
extensions of the standard model, e.g. SUSY. 

If a ferromagnetic state is formed, we would expect that the primeval plasma, where such 
bosons condensed (maybe due to a large cosmological lepton asymmetry), can be spontaneously 
magnetized. The typical size of the magnetic domains is determined by the cosmological horizon 
at the moment of the condensate evaporation. The latter takes place when the neutrino chemical 
potential, which scales as temperature in the course of cosmological cooling down, becomes smaller 
than the W mass at this temperature. However, if during the electroweak phase transition a very 
strong chaotic magnetic field was generated, the alignment of the spins of IU-bosons and the 
domain size would be determined by this magnetic field. In the course of the cosmological 
expansion such field would drop down as the cosmological scale factor squared and the spins of 
W-bosons would behave as considered above. That is, their dynamics would be determined by the 
spin-spin interactions and microscopically small magnetic domains would rearrange themselves 
into macroscopically large ones in the same way as it happens in the usual ferromagnets. 

Large scale magnetic field created by the ferromagnetism of W^-bosons might survive after the 
decay of the condensate due to the conservation of the magnetic flux in plasma with high electric 
conductivity. Such magnetic fields at macroscopically large scales may be the seeds of the observed 
larger scale galactic or intergalactic magnetic fields. This problem will be studied elsewhere, we 
only note here that the mechanism of generation of galactic or intergalactic magnetic fields is 
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unknown and presents a long standing cosmological problem, for reviews see e.g. ref. [35] • Seed 
magnetic fields generated during inflation could be quite easily uniform at galactic or intergalactic 
scales, but they are too weak, hence a huge galactic dynamo is necessary to amplify them up to the 
observed magnitude. However, it is impossible in this way to explain the existence of intergalactic 
fields. On the other hand, seed fields generated at later stages of cosmological evolution could 
be quite large (e.g. magnetic fields created at some cosmological phase transitions), but the 
characteristic scale of such fields is by far smaller than the galactic one. The mechanism suggested 
here may generate large magnetic fields and at scales which are of the order of cosmological 
horizon at the electroweak temperatures. In this sense the mechanism is unique. Still even after 
the cosmological stretching up of the characteristic wave length of the field, it remains smaller 
than the galactic radius. Nonetheless, with the "Brownian motion" reconnection of the field lines, 
the characteristic scale can be enlarged up to the galactic scale, though by an expense of their 
magnitude. Nevertheless with rather mild galactic dynamo the observed magnetic field may be 
generated. 
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